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A model for the temperature dependence of photoluminescence from 
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Photo-excited carriers, distributed among the localized states of self-assembled quantum dots, of- 
ten show very anomalous temperature dependent photoluminescence characteristics. The temper- 
ature dependence of the peak emission energy may be non-monotonic and the emission linewidth 
can get narrower with increasing temperature. This paper describes a quasi-thermodynamic model 
that naturally explains these observations. Specifically, we introduce a temperature dependent 
function to parameterize the degree of thermalization of carriers. This function allows us to 
continuously interpolate between the well-defined low and high temperature limits of the carrier 
distribution function and describe the observed anomalies in the photoluminescence spectra with 
just two fitting parameters. We show that the description is equivalent to assuming that the 
partially thermalized carriers continue to be described by equilibrium statistics, but with a higher 
effective temperature. Our treatment of the problem is computationally simpler than the usually 
employed rate equations based analyses [e.g. Sanguinetti, et al. Phys. Rev. B 60, 8276 (1999)], 
which typically also have many more under-determined fitting parameters. The model is extended 
to quantum dots with a bimodal size distribution. 
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I. INTRODUCTION 

Photoluminescence(PL) spectroscopy has perhaps 
been the most extensively used characterization tool for 
self- assembled quantum dots (QD)yJ; This is due to 
the relative ease of the measurement, the information it 
yields about the extent of quantum confinement and in- 
homogeneous broadening of the density of states and its 
very direct relevance in assessing the use of these struc- 
tures in QD lasers. 

The temperature dependent PL spectra in self- 
assembled QD have certain characteristic anomalies (Q 
d 5 S 0) At low temperatures, the linewidth of 
emission generally decreases with increasing sample tem- 
perature. The dependence of linewidth may also be non- 
monotonic with a minimum at an intermediate temper- 
ature. Secondly, the energy value corresponding to the 
peak of the emission spectrum from an ensemble of dots 
is typically observed to decrease faster with temperature 
than what is characteristic of either the bulk material or 
individual dots. It is now well recognized that both these 
anomalies are due to the non-thermal nature of the distri- 
bution of photoexcited carriers JHHEl- Energetically 
distributed bound states corresponding to QD of different 
sizes can act as trapping centers for photoexcited carri- 
ers. These electron-hole pairs, excited high up in energy 
by the pump laser can rapidly relax and get captured 
into the local potential minima in a random way. At low 
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temperatures, the finite recombination times may not be 
long enough for these carriers to cross the local poten- 
tial minima and access the lowest energy (quasi-thermal) 
state, whereas at higher temperatures the distribution 
can be expected to have better thermalized. 

Although these essential facts about the carrier 
dynamics responsible for the temperature dependent 
anomalies, are well understood l)l2j) . modeling these pro- 
cesses, even in a completely classical sense in terms of rate 
equation based models©, has turned out to be rather 
cumbersome. The rate equations l)l3j) attempt to extract 
out the most important trapping and escape processes 
with a characteristic activation energy associated with 
each© |j| HHh While this approach is both physical 
and can yield reasonable fits to the data, one typically 
requires a large number of underdetermined parameters 
(rate constants, activation energies, etc.) to explain the 
experimental data. Most of these parameters may not 
be directly accessible to experimental determination even 
in a time resolved measurement. In our opinion, the 
study of carrier dynamics is best tackle d by the more 
microscopic quantum mechanical models (|12lh while for 
describing a steady state PL experiment something sim- 
pler should suffice. 

The purpose of this paper is to attempt a physically 
motivated semi-empirical alternative that would allow for 
direct modelling of the experimental data in terms of 
small number of fitting parameters. We show that this 
is easily achieved by extending the carrier thermalization 
model ljlil) for Stokes' shift in mismatched alloy quantum 
wells. A qualitative description© Hi]) of the picture al- 
ready exists in literature. Our model is equivalent to an 
equilibrium description for localized carriers in terms of 
an effective temperature. 
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II. THEORY 

A. Formulation of the model 

Consider a semiconductor heterostructure sample, con- 
taining say, ~ 10 10 cm -2 InAs quantum dots within GaAs 
matrix, that is homogeneous on the macroscopic scale. 
The surface features appear with a size distribution, 
which when unimodal, is generally well described by a 
Gaussian functionl|T5l: Pill) . Thus the mean lateral size 
ao and the variance <r 2 /2 completely determine the mor- 
phology of a sample if, for simplicity, all the dots are 
assumed to be of similar shape with the same aspect 
ratio. Even within the assumption of a lack of phase 
coherent coupling between different quantum dots, the 
dependence of the ground state transition energy E t (a) 
on the size and shape of the dot is rendered quite non- 
trivial for such heterostructure samples. This is due to 
the finiteness of the energy barriers at materials interface, 
non-spherical shapes and the complicated strain distribu- 
tion within the dot. Furthermore since InAs is a 'narrow 
gap' semiconductor and the confinement energies can be 
even twice as much as the energy gap, one should expect 
an active participation of many bands in determining the 
energy levels. The problem of determining the size de- 
pendence of the energy levels is therefore best treated 
numerically (^3). Within a semi-empirical approach, one 
may then attempt a fit to these numerically calculated 
transition energies to a single band effective mass-like 
equation, where the confinement causes an increase in 
the effective energy gap, but now with the exact values 
of the parameters determined by a fit lfTsi) to the more 
precise numerical results. For instance, the size (basal 
length of the pyramid) dependence of the transition en- 
ergies as calculated by Grundmann et al. l)17l) can be quite 
well described by a relation of the kind 



E T {a) = E g + A/a k 



(1) 



Fixing E g at 0.42 eV, the low temperature bandgap of 
InAs, A = 2.37 eV nm 1 / 2 and k = 0.5 provided a good 
fit to the calculations in reference <|17j) . Equipped with 
a relationship between the transition energies and the 
dot sizes, one may immediately write down the expres- 
sion for the optical density of states and the absorption 
coefficient of the ensemble as a sum of absorption by indi- 
vidual quantum dots. The absorption coefficient a(E, a), 
corresponding to a single dot is essentially the product of 
the optical matrix element with the size dependent tran- 
sition energy. The density of states is therefore peaked 
at this resonance energy with a Lorentzian (19) lineshape 
due to the finite lifetime introduced by phonon and other 
scattering processes. This homogeneous broadening 7 e of 
the energies of individual quantum dots is typically less 
than a millielcctronvolt and is completely drowned by 
the ensemble inhomogeneity effects which give at least 
an order of magnitude larger contribution to the optical 
density of states. One may conveniently take the limit 
7 e — > and replace the Lorentzian by a delta function. 



The dipole matrix element \d cv \ would also be, in gen- 
eral, dependent on the effective size of the exciton rel- 
ative to the quantum dot volume. We assume ifLsj) that 
\d C v\ 2 ~ aT m . With these approximations, the ground 
state interband absorption from the complete ensemble 
quantum dots with a Gaussian size distribution of mean 
ao and variance cr 2 /2 may be written as 



= C 



m -(a-a ) /a' 



5[hv - E T (a)]da. (2) 



The constant C clubs together the normalization of the 
probability distribution function and the other size and 
energy independent parameters. 

The above integral is easily evaluated by collapsing the 
delta function after making a suitable change of variables 
and plugging in the relationship between the dot size and 
ground state transition energy given by equation ^ F° r 
the time being one may proceed without substituting the 
numerical values of k, A and E g to keep the treatment 
general. Integration of equation |2 yields (up to energy 
independent constants) 



a(E) 



(e^E s ) ( Eo -E g ) 



1/fc 



A 



Eq — E„ 



(3) 

X = (k + 1 — m)/k and Eq is the transition energy asso- 
ciated with the dot of size ao- 

In a typical low temperature non-resonant PL exper- 
iment, the carriers are generated well above the transi- 
tion energies, typically within the region of the matrix 
and then they quickly relax to the localized states within 
the quantum dots. In general, one may expect that the 
trapping efficiency of larger dots may be larger by the 
factor of a s . This should be multiplied within the inte- 
grand in equation [21 with the consequence that equation 
13 also describes the low temperature PL spectra but with 
a new value of x = (k + 1 — m + s)/k. This establishes 
the general form of the low temperature PL spectra in 
terms of the quantum dot size distribution and the corre- 
sponding exciton energies. At higher temperatures, the 
PL spectra can significantly differ due to the process of 
carrier thermalization within the localized states. Un- 
der conditions of low excitation power and in the limit of 
high temperatures, the emission spectra 1(E) is related 
to the absorption spectrum by van Roosbroeck-Shockley 
type of relationlQ), 

1(E) cx a(E,E )exp(-E/k B T), T -► oo. (4) 

ks is the Boltzmann constant and T the lattice tempera- 
ture. The proportionality constant in the above equation 
defines a temperature dependent scale factor carrying in- 
formation about non-radiative pathways. It is evident 
that the emission spectra in the low and high temper- 
ature limits differ essentially by only this exponential 
factor, that which accounts for carrier thermalization. 
Furthermore, the behavior between the two limits is ex- 
pected to be continuous and thus the essence of all car- 
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rier dynamics is contained in the extent of carrier ther- 
malization. That is, one may just rewrite the exponent 
as exp[-P(T)E/k B T] where the parameter 0(T) = 
and 1 at low and high temperatures respectively. The 
temperature dependence of the thermalization parame- 
ter, P(T) contains the essence of all the carrier dynamics 
through the multiple pathways. Within the constraint of 
< 0(T) < 1 and /3(Ti) < 0(T 2 ), if T x < T 2 the func- 
tional form of P(T) is not easy to estimate. For simplicity 
we may approximate it by a sigmoidal function 

0{T)= l + exp[-(T-7k)/Ar (5) 



T t h and A are then the only two fit parameters of the 
model. T t h is the temperature when = 0.5 and A 
corresponds to the temperature interval over which the 
transition from a 'glassy' to a quasi-equilibrium state oc- 
curs. Then, the PL spectrum at a temperature T is 



I(E,T) cx 



.4 



E-E n 



( E ~ E g) ( E 0- B s) 



IS 



I3E 



(6) 



with t — [k + 1 — m + s(l — 0)]/k. s(l — (3) is added to insure a correct interpolation between the low and high 
temperature limits, s denotes exponent of the size dependent capture cross section ~ a s of dots. (1 — (3) ensures that 
the high temperature PL is not explicitly dependent on the trapping efficiency of individual dots. The maximum of 
this equation, corresponding to the PL peak position at that temperature, is given by the solution of the following 
equation: 
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E g )i+i 



A i/k 



(7) 



where E p is the temperature dependent peak energy. 



While the values of A, E g and k were discussed earlier, 
the value of t is yet to be determined. Within the effec- 
tive mass picture, the exciton Bohr radius in (bulk) InAs 
is of the order of 35 nm. This is a bit larger than the typ- 
ical size of matured InAs/GaAs dots (height 8 nm, pyra- 
mid shape, aspect ratio=5). For spherical dots, the os- 
cillator strength \d cv \ 2 scales inversely with volume(20), 
i.e. m=3, in the limit of small dot sizes. Since our self- 
assembled dots are neither spherical nor substantially 
larger or smaller than exciton size, we do not have a 
clear understanding of the size dependence of the matrix 
element. We therefore arbitrarily assume a value of m 
that gives the simplest form to equation i.e. the value 
for which k + 1 — m — 0, with the understanding that the 
spectral shape of the absorption coefficient, being largely 
determined by the exponential function, is insensitive to 
this detail. Similarly, it is not easy to simply estimate 
the volume dependence of the exciton capture cross sec- 
tion by quantum dots and it is thus simplest to assume 
that the process of trapping of excitons by quantum dots 
it is completely random, independent of the quantum 
dot's volume. This is the usually made assumption and 
in the present context leads to a considerable simplifica- 
tion. With these assumptions |2lT) . we have put t = in 
equation Therefore Eq is the PL peak at zero temper- 
ature, corresponding to the mean sized quantum dot in 
the Gaussian ensemble. On the other hand, the shape 



of the absorption spectrum itself is no longer Gaussian 
and gets skewed toward higher energy, determined by the 
non-linear relationship between the quantum dot size and 
transition energy. 

There is also the temperature dependent contribution 
to the shift in the peak of the PL spectrum due to the 
shrinkage of the gap with temperature, which must be 
included in any analysis. Fig. [21 (dotted line) shows 
the temperature dependence of the gap (with the energy 
zero arbitrarily shifted) for InAs, with Varshni param- 
eters (a = 0.000276eFA^ 1 and 9 = 93K in equation 
below) taken from literature and therefore the position 
of the PL peak is given by 

aT 2 

E(T) = E p (T)--—. (8) 



B. Physical basis for the model 

Firstly we reiterate that the distribution function for 
non-resonantly excited carriers within an ensemble of lo- 
calized states has two distinct regimes of behavior. A low 
temperatures, the carrier distribution function is propor- 
tional to the total density of available states. At very 
high temperatures, the carriers are in thermal equilib- 
rium with the lattice. In the absence of cooperative ef- 
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fects, the crossover is expected to be smooth. The tem- 
perature dependent function (3(T), which is defined in 
Eq03 describes this crossover. 

The introduction of the parameter (3{T) amounts mak- 
ing the effective carrier temperature different from the 
lattice temperature while still assuming that the carri- 
ers are described by an equilibrium distribution func- 
tion but with a different effective temperature. The 
value of the effective carrier temperature is T/j3(T), al- 
ways greater than or equal to the lattice temperature. 
Thus the model is only approximate since it aims at an 
effective equilibrium description for an essentially non- 
equilibrium problem. In general, the idea of describing 
a non-equilibrium distribution by a higher effective tem- 
perature is not new. For example, this concept is exten- 
sively used in high electric field electron transport stud- 
ies in semiconductors. In a related study of electron lo- 
calization in Gai-adn^Asi-yNy quantum wells, we have 
recently observed that the effective temperature of im- 
perfectly thcrmalized carriers is indeed higher the lattice 
temperature|22). The concept of an effective electron 
temperature for localized carriers in disordered quantum 
wells was previously proposed in the pape r by Gurioli, 
et al. (14) and further discussed by RungeflTa) (section 
12a). Runge and coworkers fll^) have also used the equi- 
librium distribution with an effective temperature to de- 
scribe the results of their simulations of the distribution 
of partially localized excitons. Furthermore, their simu- 
lations also suggest an existence of two distinct regimes, 
which appear in our model as temperatures above and 
below T th . 

The temperature dependent behavior of the effective 
temperature T e // may, at first, appear counter-intuitive. 
The effective temperature is infinite both in the low and 
high temperature limits. This is because at low temper- 
atures, all the photoexcited carriers are frozen into the 
states they first get captured into. Thus the carrier dis- 
tribution function essentially reflects the total density of 
states. We recall that in the canonical ensemble, tem- 
perature is defined in terms of probability of occupation 
P(E) of a state of energy E through the relation, 



P{E) = 



-E/k B T 

~z ' 



(9) 



Here Z is the canonical partition function. So it follows 
that every state has an equal probability of occupation 
P(E) = l/Z when the temperature is infinite. This is 
consistent with the condition we had imposed earlier (in 
an apparently arbitrary manner) that (3 = for carriers 
to be uniformly distributed over the complete density of 
states. 

It is also tempting to compare our situation with that 
encountered in the glass transition(|23|). In context of 
glasses, the high effective temperature corresponds to an 
extra configurational entropy contribution. The notion 
of a crossover temperature is also well-established for 
glasses and is defined as the temperature at which cer- 
tain degrees of freedom get frozen. A more microscopic 



guess for the function form of the parameter (3 would re- 
quire modeling of the barriers associated with the basins 
of attraction for carriers. Very roughly, the parameters 
T e ff and A can be thought to correspond to the glass 
transition temperature and fragility respectively. 

Finally, we emphasize that the treatment presented in 
this paper is restricted to the description of the tem- 
perature dependent emission from the ground states of 
quantum dots. 



III. RESULTS AND DISCUSSION 

A. Dots with unimodal size distribution 
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FIG. 1 (a) Simulated photoluminescence spectra at different 
temperatures, llK(bottom)-160K (top) in equally spaced tem- 
perature intervals, (dots) the expected peak position assuming 
only the temperature dependent shift in the gap and (circles) 
the actual peak position. The energy difference between the 
two corresponds to the Stokes '-like shift (b) Simulated (solid 
line) and experimental (circles) PL peak energies as a function 
of temperature, (c) Simulated (solid line) and experimental 
PL full-widths at half maximum as a function of temperature. 
The experimental data was taken from reference (3) where the 
same data is fitted to an alternate model. The theoretical fit 
depicted here is with the values of a = 0.8nm, Eo = 1.22eV , 
T th = Y11K, A = 15K. 



With the above mentioned assumptions, the normal- 
ized PL spectrum is described by the following simple 
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equation (t = and k = 0.5 in equation^) 

- 1 \( A V ( A ) 

1(E) =7 (T)e ^L Vfi - E ^ Uo-« 9 ; 



(10) 



with /3(T), denned in equation^ as the fitting parameter. 

To compare with experiments, the energy dependence 
of PL intensity in equation 1 1 01 must be of course also be 
translated along the energy axis by an amount suggested 
by equation|S| 1$ (T) denotes the temperature dependent 
scale factor. While this has been previously studied and 
modelled lj24^ . we do not discuss it further because one 
usually observes the rather predictable Arrhenius-type 
activation which may as well be put in by hand. Fur- 
thermore, the non-radiative decay channels may depend 
on a particular sample's past history (e.g., the growth 
route employed to prepare the sample and the nature of 
defects) and thus not general enough to be included in a 
model like this. 

FigD shows a comparison of our calculations with the 
published results of Sanguinetti et al. (figure 4 in refer- 
ence (@)). The assumption of (3 ss at the lowest tem- 
perature fixed a — 0.8nm and E — 1.219ey. Then 
the rest of the fit was accomplished with the values of 
A = 15K and T th = YllK. This fitting procedure is 
highly constrained since it requires two separate curves 
to simultaneously fit on the basis of just two free param- 
eters. This should be compared with the rate equation 
based models, reference |&) for example had eight fitting 
parameters. 

The temperature dependent behaviour of FWHM and 
peak positions for other values of T t h and A is depicted 
in Fig|21 For ease of comparison, the ensemble character- 
istics, E and a are the same as those used to generate 
FigH an d only T t h and A are varied. The model yields 
a variety of trends at low temperatures. Some of these 
tends have not only been experimentally observed in the 
PL from self-assembled quantum dots but seem to be 
a general characteristic of emission from any inhomoge- 
neous ensemble of localized states. Thus it is tempting 
to also use the present model to analyse the PL from lo- 
calized states in the disordered quantum wells as well, in 
particular the recently much studied dilute nitrides (for a 
review see e.g., reference H^)) and InGaN (see e.g., H^)) 
quantum wells. 

But to extend this model to localized states in quan- 
tum wells would require us to guess an appropriate form 
of the density of states, which would be a combination of 
localized band-tail and higher energy extended states. In 
contrast, for self-assembled quantum dots, all states can 
be considered to be localized and the density of states can 
be more or less reliably extracted from the size distribu- 
tion which in turn is directly accessible to surface probe 
microscope analysis. Furthermore, the origin of low tem- 
perature PL in disordered quantum wells is different from 
that in quantum dots. While in quantum dots, the PL 
spectrum at zero temperature can be taken to be pro- 
portional to the total density of states (all of which are 
localized and delta- function like), only a small fraction 
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FIG. 2 Calculated temperature dependent characteristics of 
the PL spectra (a)peak positions (b)FWHM for different val- 
ues of Tth and A (in units of Kelvin) with for a fixed values 
of a — 0.8nm and Eo = 1.22eV. The dotted line in subfigure 
(a) depicts the temperature dependence of bulk InAs with the 
energy zero arbitrarily shifted. 



(corresponding to the local potential energy minima in a 
classical sense 11271) or the band tail states below the mo- 
bility edge) of the total density of states contribute to 
the low temperature PL in disordered quantum wells. 

It is also important to mention that a significant advan- 
tage of an effective equilibrium description is that one can 
continue to use the results from equilibrium theorv l)28t 
I29I) by just replacing the temperature by effective tem- 
perature T e ff. The non-monotonic behavior of the 
linewidths just follows the non-monotonicity of the effec- 
tive temperature, since the photoluminescence linewidths 
in equilibrium theory are proportional to fcgT. Sim- 
ilarly the temperature dependent "Stokes shift" lf2Sh is 
described simply by a 2 /k B T ef f ~ {k B T ef f ) 2 /k B T e f f ~ 
ksT e ff. This is evident in Fig. [21where the curves in sub- 
figures (a) and (b) seem to follow each other. Also note 
that our results are identical to the results in reference 
ffilt) in the high temperature limit of (3 — 1. 



B. Dots with bimodal size distribution 

QD ensembles often have a bimodal size distribution, 
especially during the intermediate stage of growth. This 
size bimodality also reflects in the low temperature PL 
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spectra with the appearance of two peaks corresponding 
to the two QD families. While there have been many 
experimental studies of the temperature dependence of 
the two peaks© S EJ HI S EH, the quantitative 
theoretical understanding so far has been limited to only 
the change in the relative intensity of the two peaks as a 
function of temperature ljlCt l3l|) . This was done again on 
the basis of rate equations corresponding to the change 
in the carrier population in the quantum dots belong- 
ing to the two families. Our model is easily extended to 
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FIG. 3 (Solid line) Simulated photoluminescence spectra at 
different temperatures for a bimodally size distributed quan- 
tum dots ensemble. The curves are normalized such that the 
total area of each is the same. (Dashed line) 11K spectrum 
is fitted to two Gaussians. All other curves are also well de- 
scribed as a sum of two Gaussians. Note that the higher en- 
ergy peak that virtually disappears at 220K is again visible in 
the 290K spectrum. 



describe PL in bimodally size distributed dots, by sim- 
ply accounting for the modified size distribution which 
is now a sum of two Gaussian functions, each with a 
specific mean size and size dispersion around that mean. 
Denoting by E 0l 2 , the energies corresponding to the two 
mean sizes and the (eri and a 2 ) the size dispersions in 
the two families, we may redefine the optical density of 
states which we had assumed to be proportional to the 
absorption coefficient which in turn was proportional to 



the PL spectra measured at the lowest temperature as 




i=l,2 

Here the subscripts 1 and 2 label the two families of dots 
and the ratio of the constants ao 1 2 is just the relative 
density of dots in the two families. With this new def- 
inition of ctbimod, and keeping the rest of the treatment 
the same, we have plotted simulated PL spectra at differ- 
ent temperatures in FigEI The following set of param- 
eter values was used: a± = 0.8nm, a 2 = 0.58nm, E i = 
1.22eV, E 02 = 1.28eV, T th = 150K, A = 15K. The ratio 
of the peak intensities was fixed to I(Eqi) : I(E$i) = 2:1 
at zero temperature. All these values were more or less 
arbitrarily chosen and supposed to be representative of 
the 'typical' temperature dependent behavior of PL from 
bimodally size distributed dots. 
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FIG. 4 Analysis of simulated PL spectra, Fig\^ (a) Temper- 
ature dependence of the ratio of the integrated intensity of the 
two peaks. The spectrum shows a carrier transfer to lower 
energy peak corresponding to larger dots at intermediate tem- 
peratures. In accordance with Fig\$ the peak reappears at 
higher temperatures. Note its qualitative similarity with Fig. 4 
in the work of Zhang, et al. (b) FWMH of the two peaks 

at different temperatures. 

The spectra in Fig|3|may be compared with the exper- 
imentally measured temperature dependent line shapes 
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from literature, for example Fig.l from reference 1)311) . 
Fig. 2 in reference lf32|) and Fig. 3 of reference |l"o|) . 

We note that the experimentally observed trends are 
reproduced. For ease of analysis, we again ignore the ef- 
fects of carrier loss to non-radiative decay channels and 
focus only on the effects of carrier redistribution. There- 
fore, the curves in Figdlare plotted such that the total 
integrated intensity in same in each case. The spectra 
in FigU| between 11K-150K show the transfer of carriers 
from smaller to larger dots. The intensity of the higher 
energy peak is observed to decrease. Since Tth is the sim- 
ulations was 150K, we observe that the thermalization of 
carriers is nearly complete by 200K. This leads to all the 
carriers being preferentially in larger dots and a virtual 
disappearance of the high energy peak. The carrier trans- 
fer from the smaller to larger quantum dots was found 
to be activated. Interestingly, at high enough tempera- 
tures, we observe that the high energy peak re-appears 
when the thermally broadened carrier distribution is so 
smeared that carriers have sufficient thermal energy to 
reside on smaller dots. This (somewhat unexpected) 
behavior has previously been observed in two contexts. 
Firstly, the plots of the ratio of the intensities between 
to the two peaks (Fig^ (a)) is qualitatively very similar 
to the Fig. 4 of reference flOl). Note that the analysis by 
Zhang, et al. ljlOj) was based on rate equations and em- 
ployed many more parameters to simulate a curve similar 
to FigQJa) and had no predictions about the lincwidths. 
Secondly, the reappearance of the high energy feature 
at high temperature has also been observed in context 
of pronounced appearance of the wetting layer peak|35|) 
in the high temperature PL spectra from InAs quantum 
dots. 

Finally, the temperature dependence of the positions 
and the linewidths of the two Gaussians are plotted in 
FigEJb) which shows a small but systematic decrease in 
the linewidth for the higher energy peak and a corre- 
sponding increase in the lower energy peak's linewidth. 
This feature is generic to the bimodal dots' PL. After 
150K, we observe an interesting anomaly as the behavior 
of the linewidths reverses. This 'anomaly' has also been 
experimentally observed (32) although not satisfactorily 
explained. 



IV. SUMMARY 

We have presented an extremely simple model for the 
observed anomalies in the temperature dependence of the 
PL spectra of self-assembled quantum dots. Unlike the 
previously proposed methods of analysis, we do not in- 
voke a rate equations to model the steady state carrier 
dynamics but instead introduce a physically motivated 
factor that parameterizes the degree of thermalization 
of photoexcited carriers. Since the low and high tem- 
perature limits of the carrier distribution are well un- 
derstood, the process of interpolation allows for an easy 
visualization of the carrier distribution at different tem- 



peratures in terms of this single parameter. Using this 
model, we could simulate a variety of experimentally ob- 
served trends for quantum dots ensembles with both a 
unimodal and a bimodal distribution. 
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